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Abstract 

The differential equations, which are used for matter perturbations, are usually derived from the 
Newton gravity and the Euler equation in the expanded universe. In this paper, by the explicit 
calculations of metric perturbation theory in ACDM model, we derive the general relativistic 
corrections exactly for the equations and show that the differential equations derived from the 
Newton gravity and the Euler equation are valid in the small scale region from non-relativistic 
matter dominant era onwards. However the corrections are not negligible in radiation dominant 
era. 
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I. INTRODUCTION 



In order to explain the accelerating expansion of the present universe, various models 
of the dark energy have been proposed l|-|4j]. Cosmological perturbation theory could be a 
good tool for characterizing each of the dark energy models 5, f|. Because the perturbation 
theory gives the informations about the difference from the background evolution of the 
Universe, which could depend on the characters of the models, we have a possibility to 
distinguish the characters of each model of the dark energy even if they give the identical 
evolution of the expansion in the Universe. 

The cosmological perturbation theory is often used under the sub-horizon approximation, 
which consists of the two approximations for the gravitational potential in the small scale 
a/k <C 1/H and in the Hubble scale evolution 1/dt ~ H, so that the perturbation should 
be consistent with the Newton gravity. The sub-horizon approximation is, however, merely 
an approximation so that we need to evaluate the deviations from the exact solutions in 
order to justify the approximation. While the exact calculations should be executed by 
using the general relativity, the calculation under the sub-horizon approximation gives the 
results identical with those in the Newton gravity. Recently P. Zhang has investigated the 
general relativistic corrections from the Newton gravity and has shown the corrections could 

n 

not be neglected in the large scale region [7(. As we show in Section IHIt the corrections 
are negligible in the small scale region by the exact evaluation of the general relativistic 
corrections from matter dominant era onwards. But we also show that the corrections are 
not negligible in radiation dominant era even if we consider the matter perturbations in the 
small scale region. 

We use units of &b = c = K — 1 and denote the gravitational constant 8nG by n 2 in the 
following. 



II. COSMOLOGICAL PERTURBATIONS 



In this section, we consider the perturbed equations in ACDM model by using the New- 
tonian gauge. We assume, however, that the matter components are given by a single 
barotropic fluid. Therefore the following arguments can be applied only in radiation domi- 
nant era and from non-relativistic matter dominant era onwards. 
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We begin with the Einstein equation including the cosmo logical constant: 

Rtw - -^QiivR = -#^A + k 2 T^ v . (l) 

If we assume that the spacial curvature should vanish, the Friedmann-Lemaitre-Robertson- 
Walker (FLRW) equations are given by substituting the FLRW metric ds 2 = —dt 2 + 
5 ij a(t) 2 dx i dx j into Eq. ©: 

3H 2 = k 2 P + A, 

-a% (2H + 3H 2 ) = K 2 a 2 6 ijP - a 2 %A, (2) 



where we assume the energy momentum tensor of perfect fluid, = pg^u + ip+p) 



0, uq = —1, where is the four vector satisfying g^u^Uy = —1, p is the energy density, and 
p is the pressure. By substituting the FLRW metric with perturbation in the Newtonian 
gauge, 

ds 2 = (-1 + 2<$>)dt 2 + 5 ija (t) 2 (l + 2^)dx i dx j , (3) 
into the perturbed Einstein equation, 



we obtain the following equations in the leading order of the perturbation, 

k 2 

-GH 2 $ - 2—^ - 6Hd ^ = - K 2 5p, (5) 
or 

2di(H$ + doV) =K 2 (p + p)6 UlJ (6) 
a" 2 ^($-vl>)=0, (i^j), (7) 



2H8 - AH - 6H 2 ) $ - — + -Y + 2<9 a + 6H8 * 

=K 2 8p, (not summed with respect to i) , (8) 

where Eqs. (JSJ, (EJ), ([7]), and (jSJ) express the (00), (Oi), (ij) for i ^ j, and (ij) for « = j 
components of the Einstein equations, respectively. In fl5]) and k represents the wave 
number which appears from the derivative with respect to the spacial coordinates (k 2 = 
—djdj) by the Fourier transformation. We treat the energy momentum tensor as that of the 



3 



perfect fluid so that the perturbations of the energy momentum tensor are given by 

5T° = -8p, (9) 

ST° l = (p + p)Su i , (10) 

6T = -a- 2 (p + p)5u u (11) 

6T) = SijSp, (12) 

where we use 5uq = 8 goo/2, which is obtained from the condition g^u^Uy — — 1. When we 
consider the scalar perturbation, we decompose 5ui as di5u + 8uJ , where d % 5u( = and 
put 8u( = and keep only the scalar part 8u. Then 5ui is expressed as 5ui = di5u. The 
equation of the perturbation for the matter is given by considering the fluctuation in the 
equation of continuity V M T^ = 0, 

5p + 3H{5p + Sp) + a- 2 cU(p + p)5u i ] + 3^(p +p) = 0, (13) 
a~ 3 d {a 3 (p + p)5ui\ + - (p + p)^$ = 0. (14) 

Basically, we obtain the time evolution of matter density perturbation 5 = 8p/p from the 
above these equations. It should be noted, however, there are unknown variables 5ui and \& 
and therefore we must need to use the other equations such as (jSJ), ©, and so on. 



III. SUB-HORIZON APPROXIMATION AND THE EXACT EQUATION 
A. Deduction of the differential equation 

In the following, we choose \l/ = $ in Eq. (JTJ). If we use the sub-horizon approximation, 
which is given by H 2 ^, Hd ^, dfa < k 2 ^/a 2 , Eq. (jSJ) gives, 

k 2 

2—m ~ K 2 5p. (15) 
a 1 

This reproduces the Newton potential if we replace k 2 /a 2 and 8p by —(f), —V 2 and p, 
respectively. Such an approximation is correct when we consider the case of small sound 
speed and the region that the expansion speed of the Universe is negligible. But we should 
consider the calculation error without any assumptions. Thus we now evaluate the correction 
by the general relativity to the differential equations of matter perturbation. When we use 



4 



the sub-horizon approximation and when w 
of matter perturbation is given by, 



p = dp = 0, we find the differential equation 



6 + 2H5 - -fl m H 2 6 ~ 0, 



(16) 



which is given by using Eqs. ([7j), ffTB"]) . (TH|) . and f|T5|) . Here we define 5 = <9o<5, 5 = <5p/p, 
and Q m = n 2 p/3H 2 . The equation (flBl) is not changed from the equation given by using the 
Newton gravity and the Euler equation [8|. If we represent Eq. (TT6]) by using dimensionless 
variable N = In a, we obtain 



d 2 5 
dN 2 



1 



d5 



o - o^'cfr ] - ^ 0, 



(17) 



where w c s = — 1 — 2H / (3H 2 ) is defined by Eqs. in (j2J). On the other hand, we use Eqs. (JSj), 
([6]), ([7]), and ffTBl to obtain the exact equation of matter density perturbation. First, Eqs. 
(JSJ) and (|7j) give the one dimensional differential equation for 



H 



Furthermore, this equation (TT5]) and Eqs. ( TT3"]) . 



5p 



3£T(1 + c£) - (1 + «;) 
I- <| 3(1 + w)Hp + (1 



3aH 
, and (j7|) give 
K 2 p fc 2 



(18) 



2# 3a?H 
, pfc 2 2A; 4 



5p 



a 2 H 3a i Hn 2 ' 

Here the sound velocity c 2 is defined by 5p = c 2 5p. By differentiating Eq. ( |19p with respect 
to N and eliminating the terms proportional to \I/ and \l/ by using ( fl8l) and ( TT9l) . we find the 
differential equation expressed only by the variable Sp. If we use the dimensionless variables 
5 and N = In a, the differential equation is given by 



rf 2 £ 



+ 1 



1 + w)tt m + 3(c; 



hi 



_d_ 

dN 

k 2 
3a 2 H 2 

k 2 



2k 4 



3(1 



(2 



3a 4 H 2 n 2 p 
- 3w — 3c 2 + 3w e ff) 



3a 2 H 2 



3(w — c 



;(1 




0. (20) 
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This is the exact differential equation of matter perturbation in the Newtonian gauge. If we 
use only small scale approximation, H 2 <C k 2 /a 2 in order to compare Eq. (120]) with Eq. ffTTj) . 
Eq. (120 p has the following form: 



d 2 6 
dN 2 



3w cS + 3w — -(1 + w)VL rn + O 



k 2 



+ 



wk 2 



3 -a 

2 V 



w)(l + 3w)Q m + O 



a 2 H 2 

k 2 
a 2 H 2 



-i 



d6_ 

dN 

5 = 0. 



(21) 



^a 2 H 2 

Here we set c 2 = w for simplicity. To derive the above equation f[2"Tj) . we used the following 
expansion, 

k 2 



d_ 

dN 



In 



2k 4 



3a 4 H 2 K 2 p 
3w c r + 3w -+ 



3 1 



, 3a 2 n 2 p ^ 
'l + w )(l + 3w)^f+0 



k' 2 



a 2 H 2 



(22) 



We need to take w = to compare with Eq. ( Tl71) because Eq. ( TlT|) can be only applied for 
the non-relativistic matter dominant era onwards. 



d 2 5 
dN 2 



2 " 2 Weff 



A- 2 



o 



a 2 H 2 



-i 



5 = 0, (23) 



where we use the equation (1 + w)VL m = 1 + w c q given by (jSJ). Equation f l23|) is a same 
equation with Eq. flT7]) up to O ({A; 2 /(a 2 i7 2 )} ). So that it is shown that the Eq. flUD is 
valid in the small scale region H 2 <C k 2 /a 2 . 

Whereas the scale we are interested in isn't always deep inside the horizon in radiation 
dominant era. Therefore we need to consider the whole corrections in Eq. ( 1201) . Even 
if we consider the region deep inside the horizon, we can see that the general relativistic 
corrections are not negligible in radiation dominant era from Eq. (|2T|) . If we express Eq. 
( 12T|) by using time derivative, we obtain 

k 2 x-i N 



2 - 3w + O 
wk 2 3 



a 2 H 2 



HS 



a?H 2 



w ){i + 3w)n m + o{^ f2 



5 = 0. 



(24) 



This equation (|24|) is equivalent to the equation in 0, [uj] which is derived by the different 
ways in general relativity. 
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B. Consistency check 

We obtain the differential equation of ^ from Eqs. (Tl8|) and ( TT9|) in a same way as we 
obtain the differential equation of 5 in the last subsection. 

+ I 1 + 3 ( c * - + + w)Q m j ^+{^-2+ 3(c 2 - u,^) j * = 0. (25) 

The differential equation for \1/ can be, however, also obtained by another way Eliminating 
the terms proportional to 5 and $ by using Eqs. (jSJ), (EJ), and (jHJ), we obtain the following 
equation, 

d 2 * /5 2 3 \ dm f c 2 *; 2 . 2 ,1 T 

^ + (5 + 3cJ - ^ + + 3(c* - ^)) * = 0. (26) 

The obtained expression of Eq. ( 126]) might appear to be different from that in ( 1251) but if 
we use the equation (1 + w)Q m = 1 + w c g, we can find Eq. f f26|) is surely identical with 
(|25|) . This might be considered to be a trivial result but in case that there are non-trivial 
contribution from the diagonal parts of the metric perturbation, the obtained results are 
not identical with each other. 



IV. CONCLUSION 

In this paper, we derived the exact differential equation of matter perturbation in the 
Newtonian gauge. This equation ([211]) gives the standard equation of matter density pertur- 
bation, 5 + 2H5 — 3Q m H 2 5/2 ~ 0, from non-relativistic dominant era onwards by the small 
scale approximation H 2 <C k 2 la 2 . Whereas in the radiation dominant era, Eq. ( 1201) gives the 

fin 

same result in [9|, [lfj at the deep inside the horizon. So that the Eq. (120]) can be thought as 
the generalization of those equations. In particular Eq. (I20p is useful in radiation dominant 
era because the small scale approximation could not work well in radiation dominant era. 
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